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KeksinnOn kohteena on menetelma puustotunnusten maarittamiseksi puuston 
ylapuolella mittaavan laitteen avulla. Menetelmassa puustosta kerfitaan kolmiulot- 
teista tietoa puustoalueelta kayttamaiia sellaista maarSa naytepisteita, etta yksit- 
taiset puut tai puuryhmat erottuvat toisistaan. KerStyista tiedoista tuotetaan 
kolmiulotteinen puuston pituusmalli, josta maaritetaan puustotunnuksia, jotka ovat 
yksittaisten puiden tai puuryhmien ominaisuuksia ja/tai nfiden avulla saatuja laa- 
jempien alueiden tietoja. Keksinn6n kohteena on my6s tieokoneohjelma, jolla 
voidaan toteuttaa menetelman toinen ja kolmas vaihe. 




4, 



Referansftl- 
uama 



112402 

2 

1991. Satellite image-based national forest inventory of Finland. International 
Archives of Photogrammetry and Remote Sensing. 28: 419-424). Taliein tyypillisesti 
maastoinventoinnin koealatuloksia kaytetaan menetelman opetusaineistona ja koko 
satelliittikuvaa kaytetaan yleistamaan nama tarkasti koealoilta mitatut tiedot koko 
5 kuvan alueelle. Edellytyksena onnistuneelle ratkaisulle pienilla alueilla on se, etta 
satelliittikuvasta tulee loytad piirteitd (kanavia, niiden suhteita tms), jotka korreloivat 
voimakkaasti koealoilta maSritettyjen puustotunnusten kanssa. Niinpa menetelma 
onkin sovelias vain erittain laajojen alueiden inventointiin. Parantamalla 
kaukohavainnointiaineistoa, esimerkiksi lentokoneesta tehtavilla kuvauksilla, voidaan 

10 parantaa myos sen tarkkuutta. Huolimatta tastS kuvioittaisen arvioinnin vaatimaa 
tarkkuutta (noin 15% virhe sallittu, R. PSivinen, A. Pussinen, and E. Tomppo, 1993, 
"Assessment of boreal forest stands using field assessment and remote sensing", 
Proceedings of Earsel 1993 Conference "Operationalization of Remote Sensing", 
ITC Enshedene, The Netherlands, 19-23 huhtikuu, 1993, 8s.) ei ole 

15 kaukohavainnointimenetelmin saavutettu. 

Esimerkiksi Suomessa metsaorganisaatioiden toimesta tehtava kuvioittainen 
inventointi tehdSSn kokonaan maastotyona ja valtakunnan metsien inventointiin 
kayteta§n satellittikuvia (kuten Landsat TM-kuvia, joiden erotuskyky on 30 m). Varsin 

20 perusteellinen selvitys eri aineistojen tarkkuudesta kuvioittaiseen arviointiin on 
esitetty julkaisuissa {J. HyyppS, HyyppS, H., Inkinen, M., Engdahl, M., Linko, S. and 
Zhu, Y-H., 1998. Accuracy of different remote sensing data sources in the retrieval of 
forest stand attributes. Proceedings of the First International Conference on 
Geospatial Information in Agriculture and Forestry. Lake Buena Vista, Florida, USA, 

25 1-3 kesSkuu 1998, Volume I, s. 370-377. ja J. Hyyppa, Hyyppa, H., Inkinen, M, 
Engdahl, M., Linko, S., Zhu, Y-H., 1999a, Accuracy comparison of various remote 
sensing data sources in the retrieval of forest stand attributes, Journal of Forest 
Ecology and Management (in press)). Hakijan tiedossa ei ole olemassa aiemmin 
tunnettuja kuvioittaisen arvioinnin vaatimuksia tayttavaa kaukohavainnointiaineistoa. 

30 

Toisenlainen menetelma metsapalstan runkosuureiden laskemiseksi on esitetty Fl- 
patentissa 101016 B. Menetelma optisesti rekisteroi valitussa pisteessa tietylla 
sateelia kaikki puut 360 asteen kierrolla. Mittauksessa kaytetaan esimerkiksi 



3 



112402 



AccuRange 3000-LIR Laser etaisyysmittaria, jonka mittaustarkkuus on 1/65535-osaa 
360° :sta tai pulssiahturia. Rekisterointi suoritetaan kiertokulman mittarilla siten, etta. 
tangenttikulma lasketaan absoluuttisten kiertokulmien erotuksena rungon molemmin 
puolin havaituista etaisyysmittauksen epajatkuvuuskohdista. Menetelmaa (Fl-patentti 

5 101016 B) voidaan kayttaa perinteisen koealatietojen automatisointiin, mutta se 
edellyttaa silti tybskentelya maastossa ja on varsin hidas (yhden kierroksen 
mittausajan tulee olla 1-6 min, jotta tietokone ehtii siirtaa mittaustiedot). Menetelma. 
perustuu myos lapimittojen kautta tapahtuvaan inventointiin, joka tapahtuu 
horisontaalisella mittauksella. Ongelmana tassS Fl-patentin 101016 B mukaisessa 

10 ratkaisussa on etenkin se, etta esitetty menetelma hitautensa puolesta soveltuu vain 
pieneen nayteenottoon koko puustosta. 

Lentokoneesta ja helikopterista on mitattu metsaaluiden puuston pituutta 
aikaisemmin myos kayttaen lasereita ja radiotaajuustutkia. Aikaisemmat tutkimukset 

15 perustuivat poikkileikkauksen mittaamiseen metsa.alueilta (lentolinjaa pitkin mitattiin 
puuston pituutta laserin/radiotaajuustutkan keilan valaisemalta alueelta). 
Esimerkkeja tailaisista tutkimuksista ovat mm. R. Nelson, Krabill, W.B., and Maclean, 
G.A. 1984, "Determining forest canopy characteristics using airborne laser data", 
Remote Sensing Environment, 15:201-212, ja_ J. Hyyppa, Hallikainen, M., 1996. 

20 Applicability of airborne profiling radar to forest inventory. Remote Sensing 
Environment, 57: 39-57. Yksittaisia puita ei ole naissa tutkimuksissa pystytty 
analysoimaan, koska kuvat olivat kaksiulotteisia poikkileikkauksia. Mittauksilla on 
saatu puuston pituus ja siita regressiokaavoilla johdettua muita pituuden 
johdannaisia. Mm. tilavuuden osalta puusto tilavuuden virhe oli parhaimmillaan 26.5 

25 %, mik§ ei riita operatiiviseen mittaukseen. 

N§sset (esim. E. NSsset, "Determination of mean tree height of forest stands using 
airborne laser scanner data", ISPR J. Photogramm. Remote Sensing, 52, s. 49-56, 
1997.) on pystynyt tuottamaan laserkeilauksella tasaisesti naytteita metsasta, mutta 
30 puustotunnusten estimointi tapahtui puhtaasti tilastollisesti likimain samalla tavoin 
kuin aikaisemmat profiilimittaukset. Mm. keskipituusestimaattina kaytettiin 
tietynkokoisen ikkunan sisalta otettujen lasernaytteiden minimeja ja maksimeja. 
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Vuonna 1999 Hyyppa ym. (J. Hyyppa, Hyyppa, H., Samberg, A., 1999, Assessing 
Forest Stand Attributes by Laser Scanner, Laser Radar Technology and Applications 
IV, Proceedings of SPIE, 3707: 57-69.) esitti, etta on mahdollista mitata tarpeeksi 
tihe§lla pulssiverkolla metsan valtapuuston pituus. Tassa tutkimuksessa kokeiltiin 
5 myos pituusnaytteiden avulla tapahtuvaa tilavuuden estimointia samalla tavoin kuin 
edella mainituissa pituusmittauksissa, joissa kaytettiin laserlaitteita ja 
radiotaajuustutkia. LisSksi tutkimuksessa esitettiin todellisuutta vastaava 
laserkeilaindatan avulla tuotettu puuston pituusmalli. Tutkimuksessa ei tunnistettu 
yksittaisia puita eika puustotunnuksia. 

10 

Myos Gunilla Borgefors ym. (Gunilla Borgefors, Tomas Brandberg, Fredrik Walter 
"Forest parameter extraction from airborne sensors", APRS, Vol.32, Part3- 
2WS,"Automatic Extraction of GIS Objects from- Digital Immaging", Munchen 8-10 
syyskuuta 1999, s. 151-158) on esittanyt laserdatan hyvaksikayttamista puiden 
15 tunnusten mSarittamisessa. Julkaisussa on maaritelty puiden lukumaara ja niiden 
latvojen koko lasertiedon avulla saadusta kuvasta analysoimalla kuvasta nakyvia eri 
alueiden korkeuseroja referenssitasoon nahden. Yksittaiden puiden todellisia 
pituuksia ei tassa julkaisussa analysoitu. 

20 Laser- ja tutkatekniikkaa on myos yleisesti kSytetty maastomallien tekemiseen 
mittaamalla ylhSSItapdin lentokoneesta kohteen ja tutkan (yleisnimitys seka 
lasertutkalle etta radiotaajuustutkalle) vSlinen etaisyys pulssin kulkuajan perusteella. 
Lasertutkaa kutsutaan yleisesti myos lidariksi (light detection and ranging). Talloin 
laserkeilain pyyhkaisee tutkan muodostamaa laser-pulssia lentosuuntaa vastaan 

25 kohtisuoraan. TSIIS tavalla katetaan koko kohteen ala. Alueelta saadaan lahes 
vierekk^isia keiloja (keilan halkaisija tyypillisesti muutama kymmen senttimetria), 
joista jokaiselle saadaan maaritettya x-, y- ja z-koordinaatit. Naita pisteita 
analysoimalla voidaan tuottaa erilaisia maastomalleja. Radiotaajuustutka toimii 
samalla periaatteella kuin lasertutka, lahetettavan signaalin taajuus on vain 

30 radiotaajuusalueella. Radiotaajuustutkilla keilaus on mahdollista toteuttaa joko 
elektronisesti tai mekaanisesti. Tyypillisesti radiotaajuusalueella antennin 
synnyttaman keilan koko maanpinnalla on useita metreja, toisaalta SAR-tekniikalla 
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(engl: Synthetic Aperture Radar) voidaan parantaa huomattavasti radiotaajuustutkan 
erotuskykya. Radiotaajuustutkilla voidaan mitata samanaikaisesti seka latvan etta 
maanpinnan etaisyys tutkaan nahden, joten erilaisten maastomallien tuottaminen 
radiotaajuustutkilla on yksinkertaisempaa, toisaalta erotuskyky on vaikea saada yhta 
5 hyvaksi kuin nykyisilla ja tulevilla laserkeilaimilla, jotka pystyvat lahes 100 kHz 
pulssintoistataajuuteen, jolloin kohde voidaan kartoittaa jopa 0,5 metrin 
erotuskyvylla. 

Keksinndn tarkoituksena on aikaansaada kaukohavainnointimenetelma metsien 
10 inventointiin ja puustotunnusten maarittamiseen, joka on aikaisempia 
kaukohavainnointimenetelmia nopeampi, tarkempi, mahdollistaa paremman tulosten 
jatkokasittelyn ja on suhteellisen kustannustehokas. 

15 YHTEENVETO KEKSINNOSTA 

Keksinndn mukainen menetelma puustotunnusten maarittamiseksi puuston 
ylapuolella mittaavan laitteen avulla on paaasiassa tunnettu siita, etta puustosta 
kerataan kolmiulotteista tietoa puustoalueelta kayttamalla sellaista maaraa, 
20 naytepisteita, etta yksittaiset puut tai puuryhmat erottuvat toisistaan. Keratyista 
tieddista tuotetaan kolmiulotteinen puuston pituusmalli. Saadusta pituusmallista 
maaritetaan puustotunnuksia, jotka ovat yksittaisten puiden tai puuryhmien 
ominaisuuksia ja/tai niiden avulla saatuja laajempien alueiden tietoja. 

25 Puustotunnusten maarittaminen keksinndssa tapahtuu edullisimmin kolmessa 
vaiheessa 

1) Hyvaresoluutioisen kolmiulotteisen mittauksen tekeminen puustoalueelta. 

2) Kolmiulotteisen puuston pituusmallin laskeminen tuotetusta aineistosta. 

3) Puustotunnusten maarittaminen puuston pituusmallista yksittaisten puiden ja 
30 puuryppaiden avulla. 

Keksinndn mukaisen menetelman toinen ja kolmas vaihe voidaan toteuttaa 
tietokoneohjelman avulla. 
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20 



Hyvaresoluutioisella tai hyvan erotuskyvyn omaavalla tassa tarkoitetaan aineistoa, 
josta pystyy havaitsemaan erikseen yksittaisia puita. 

Pohjoisella havumetsavyohykkeella, monilla talousmetsaalueilla ja muissakin 
metsaalueilla on selvia aukkoja puiden valissa. Esimerkiksi tiheassa metsassa . 
Suomessa yli 30 % laserpulsseista kuitenkin heijastuu maasta. Lisaamalla 
laserkeilaimien pulssitaajuutta (lahetettyjen laserpulssien maaraa sekunnissa) 
merkittavasti voidaan laserkeilaimen avulla saada naytteita jokaisen puun latvasta ja 
puiden valisesta maasta. Laserpulsseja tulee talloin olla useita neliometria kohden. 
Tama mahdollistaa hyvaerotuskykyisen puustokartan tekemisen laserkeilaimen 
tuottamasta aineistosta. Kun tata aineistoa kasitellaan, siita saadaan laskettua 
maastomalli ja latvamalli seka niiden erotus, eli puuston pituusmalli. Pituusmallia 
analysoimalla esimerkiksi hahmontunnistuksen menetelmin voidaan maarittaa 
yksittaisten puiden sijainti, puiden pituus, latvuksen leveys, puulaji ja niiden avulla 
laskea kullekin puulle lapimitta, runkoluku, ika, kehitysluokka, pohjapinta- ala seka 
tilavuus. Koealoille ja kuvioille voidaan vastaavasti laskea kyseiset tiedot 
puulajiluokittain. Keksinnossa on myQs esitetty, miten vanhaa puustotietoa ja 
tietamysta (tietamyspohjaisia jarjestelmia) voidaan kayttaa parantamaan 
estimoitavien puustotunnusten tarkkuutta. 



TSssa keksinnossa esitetty menetelma on tarkin kaikista tahanastisista metsien 
inventoinnin kaukohavainnointimenetelmista. Keksinnon etuna tyypillisiin 
kaukohavainnointimetelmiin nahden on my6s se, etta keksinnon menetelmalla 
mitataan kohteesta selvia fysikaalisia tunnuksia etaisyysdatan muodossa, joille 
5 voidaan johtaa yleisesti tunnettujen yhtaioiden avulla puun tunnuksia, kuten tilavuus. 
Nain ollen menetelma ei valttamatta vaadi koealojen kayttda opetusaineistona, mika 
pienentaa menetelmasta johtuvia kustannuksia. 

Hyvaresoluutioisen kolmiulotteisen kuvan tuottaminen on edellytys menetelman 
0 kaytolle. Mittauksien tiheys (perakkaisten ja vierekkaisten pulssien/naytteiden 
etaisyys maanpinnalla) vaikuttaa voimakkaasti siihen, mita laserin tuottamasta 
aineistosta on mahdollista maarittaa. . Keksinnossa on oivallettu hyodyntaa 
lasertutkien kasvanut pulssitoistotaajuus menetelmassa, joka perustuu enemman 
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hahmontunnistukseen (saadaan useita osumia kohteeseen neliometria kohden, 
jolloin useiden kohteiden piirteet alkavat nakya) kuin nykyisen tekniikan tason 
mukaiseen puhtaasti tilastolliseen laskentaan. Puissimaaran kasvattaminen 
keksinnon mukaisen menetelman tavoin tuo merkittavaa lisSetua metsien 
5 inventointiin. Suuri nSytteiden maara mahdollistaa hyvan erotuskyvyn ja talloin on 
mahdollista mSaritt§3 yksittaisen puun ominaisuuksia. Aiemmin laserdatan avulla ei 
ole tuotettu kuvia, joista on pystytty erottamaan ja johtamaan tarkasti yksittaisten 
puiden todelliset ominaisuudet, joista mm. tilavuus voidaan laskea. 

10 Keksinnon mukainen menetelma pystyy yleisesti ottaen tuottamaan puuston 
tilavuusestimaatteja pienemmalla kuin 15 % virheella kuviotasolla pohjoisella 
havumetsavyfihykkeella. Menetelma on sovellettavissa myos muuallakin, erityisesti 
tropiikin istutusmetsissa. Tama on mahdollistuu siten, etta on mahdollista tehda 
ylhaalta nakyvasta puustosta puustokartta ja koska kartta on suoraan digitaalisessa 

15 muodossa, ja nykyaikana suurin osa metsatiedoista on paikkatietojarjestelmissa, se 
mahdollistaa jopa yksittain puin alueen seuraamisen ja toimenpiteiden paivittamisen, 
mikS esimerkiksi puisto- ja arvoalueilla voi olla tarpeen. Menetelman 
kustannustehokkuutta nostaa jatkuva laserkeilaimien ja tutkatekniikan kehittyminen. 
MenetelmSn automatisointi (keratyn aineiston prosessointi tietokoneella) on myos 

20 mahdollista. 

Seuraavassa keksinto esitetaan yksityiskohtaisesti kuvioiden ja esimerkkien avulla, 
joita ei ole tarkoitettu keksinnon rajoittamiseksi millaan tavalla. 

25 

KUVIOT 

Kuvio 1 on yleiskuva keksinnon mukaisen menetelman mittaustilanteesta. 
Kuvio 2 esittaa laserkeilaimen avulla saatavaa naytteenottotiheytta metsaisesta 
30 kohteesta. 

Kuvio 3 esittaa naytteiden maaran ja erotuskyvyn valista yhteytta metsaisella 
alueella. 
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Kuvio 4 on todellisella laserkeilaindatalla laskettu minimipinta metsaalueelle (kunkin 
pikselin sisalle kuuluvista pisteista on valittu z-arvoltaan pienin). Pikselikoko 
50 cm, koko alue 125 m x 125 m. Puuttuvat arvot koodattu 150:ksi. 

Kuvio 5 on todellisella laserkeilandatalla laskettu maksimipinta metsaalueelle (kunkin 
pikselin sisalle kuuluvista pisteista on valitty z-arvoltaan suurin). Puuttuvat 
arvot koodattu 0:ksi. 

Kuvio 6 on kuvioiden 4 ja 5 datalle laskettu maastomalli. Maastomalli kuvaa suoraan 
maaston korkeutta kyseisessa pisteessa meren pintaan nahden. 

Kuvio 7 on kuvioiden 4 ja 5 datalla laskettu puuston pituusmalli. Pituus pituusmallista 
pystytaan todellakin erottamaan yksittaisia puita ja puurykelmia. Kuvan 
erotuskyky on 0.5 m. 

Kuvio 8 on esimerkki segmentoiduista latvuksista. Segmentoinnilla maaritetaan 
latvuston reunat. 

Kuvio 9 on esimerkki yksittaisten puiden pituuden mittaustarkkuudesta keksinnon 
osoittamalla tavalla. 

Kuviot 1-3 ovat karkeita kuvioita mittaustapahtuman selkeyttamiseksi. Kuvio 8 on 
karkea kuvio latvusten aariviivojen maarittamisesta. 



KEKSINNON YKSITYISKOHTAINEN KUVAUS 

Hyvaresoluutioisen kolmiulotteisen mittauksen tekeminen puustoalueelta 

Keksinnossa tutkalaite, jota kaytetaan mittauksissa, voi olla joko levean keilan 
muodostama laserkeilain tai radiotaajuinen keilaava tutka, jotka ovat toiminnaltaan 
samanlaisia. Tarkeinta on tuottaa tutkaperiaatteella hyvaresoluutioista 
kolmiuloitteista tietoa metsasta, josta yksi puin tai puuryhmin puustotunnisten 
maarittaminen tehd§an. 

Mikali laserkeilaimen synnyttaman keilan leveys on luokkaa 1 m, laite voi 
rekisteroida saman pulssin sisalla puuston aiheuttaman kaiun. Samasta kaiusta voi 
talloin maarittaa seka maanpinnan ja latvan etaisyyden tutkasta. Keilaavalla 



112402 



radiotaajuustutkalla tilanne on samanlainen. Ero kapean keilan muodostavaan 
laserkeilaimeen on talldin vain, se, etta myds maakaiulle saadaan useammin 
mittaustieto, eli signaalin tunkeutumiskyky puuston lavitse on parempi. Tallaisilla 
laitteilla on siten helpompi maarittaa luotettavasti maastomalli. Kapean keilan 

S muodostava laserkeilain edustaa talla hetkella olevaa potentiaalista instrumenttia 
puuston pituusmallin msaritt§miseen, joten jatkossa keskitytaan kuvaamaan, miten 
silla saadaan puuston pituusmalli tuotettua ja miten sen aineistoa voidaan hyddyntaa 
metsien inventoinnissa. Vaikka jatkossa siis puhutaan laserkeilaimesta (kapean 
keilan muodostavasta), keksintd kattaa myds yleisella tasolla periaatteet, miten 

10 levean keilan muodostavat laserkeilaimet ja radiotaajuiset keilaavat tutkat toimivat 
metsaolosuhteisa. Jalkimmaisilla on vain mahdollisuus tuottaa raakatuloksena 
kolmiulotteista tilavuusmallia kuin edellinen tuottaa kolmiulotteista pintamallia. 

Kuviossa 1 on esitetty yleiskuva keksinndn mukaisen menetelman mittaustilanteesta. 

15 Mittaus tehdaan lentavasta aluksesta 1 (lentokone, helikopteri, miehittamatdn 
lentoalus, lennokki) puuston ylapuolelta. Naista lentokone on talla hetkella tasaisilla 
alueiden paras alusta, helikopteri vuoristoisilla alueilla ja kustannuksiltaan edullisin 
on automaattiohjattu miehittamatdn lentoalus tai lennokki. Laseriaite 2 koostuu 
keilainosasta, joka aiheuttaa lentosuuntaa vastaisen poikkeutuksen, lasertykista, 

20 joka tuottaa laserpulssit 3, ja ilmaisinosasta, joka tulkitsee vastaanotetun signaalin ja 
maarittaa sen perusteella etaisyyden kohteeseen. Kuvion 1 lentokoneessa 1 olevan 
laser-keilaimen lasertykin lahettama laserpulssi 3 osuu kohteeseen ja palaa takaisin 
ilmaisimelle, jolloin voidaan maarittaa kohteen ja laserin valinen etaisyys pulssin 
kulkuajan perusteella. Kun lasertykin asento ja paikka ovat tarkasti tunnetut, mitattu 

25 etaisyys lasertykin ja kohteen 4 tassa tapuksessa puun) valilla voidaan muuttaa 
korkeudeksi; tahan perustuu koko laserin avulla tapahtuvan pintamallin mittaaminen. 
Laserlaitteen 2 asento ja sijainti maaritetaan tyypillisesti inertiajarjetelman ja GPS- 
mittauksen avulla (Global Positioning System, maailmanlaajuinen satelliittiverkosto 
paikantamiseen). Inertiajarjestelma mittaa joko pelkkaa asentoa tai asentoa ja 

30 paikkaa inertia-antureita kayttaen. Laserkeilain tekniikan yhteydessa puhutaan 
yleensa lyhenteesta INS, joka on taydellinen 6-ulotteinen mittausj§rjestelma (seka 
asento ja paikka). Laserkeilaintekniikan GPS-mittauksessa on yksi GPS-vastaanotin 
5, joka on samassa lentavassa alustassa kuin laserkeilain ja maassa on toinen GPS- 
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vastaanotin referenssiasemana 6 tutkittavan alueen lahietaisyydessa, esimerkiksi 30 
km sateelld. 

Koska laser ja radiotaajuustutkat ovat aktiivisia mittalaitteita (Iahett§vat signaalia, 
5 jonka vastaanotetaan kohteesta heijastuttuaan), voidaan mittauksia tehda myos 
yollS. 

Kuvan tuottaminen tapahtuu keilaimen avulla, joka pyyhkaisee laser-pulssia 
lentosuuntaa vastaan kohtisuoraan. Toisessa suunnassa kuva muodostuu lentavan 

10 alustan siirtyessa pitkin lentolinjaa. Talla tavalla katetaan koko kohteen ala. 
Lasertykin ja keilaimen ominaisuuksista (pulssitoistotaajuudesta, rinnakkaisten 
pulssien maar§st§, keilaimen maksimaalisesta keilauskulmasta ja keilaustavasta), 
lentokorkeudesta ja lentonopeudesta riippuen maaraytyvSt mm. maanpintaan 
vierekkaisten keilojen (keilan halkaisija tyypillisesti muutama kymmen senttimetria) 

IS etaisyydet. Keilausmekanismeja on useita. Kartiokeilauksessa mittauskulma 
kohteeseen pysyy koko ajan vakiona. Nk. pushbroom-keilauksessa rinnakkaisten 
keilojen orientaation toisiinsa nahden on aina vakio. 

Kuviossa 2 on esitetty tuloskuvio siitd, miten yksittaisilla keiloilla 3 peitetaan 
70 mitattava alue mittauspisteiden 7 avulla. Mitattujen pisteiden 7 sijainnit muodostavat 
kuviossa 2 saannollisen kuvion lentotieta pitkin. Viitenumero 8 esittaa puiden 
latvuksen peittaman alueen. 

Metsien mittaamisen kannalta tarkeata on tehda mittaukset mahdollisimman 
25 kohtisuoraan, jotta katvealueita ei syntyisi puiden taakse. Mittauskulma tulee olla 
mielellaan pienempi kuin 10 astetta vertikaalista. Koemittauksissa on havaittu 
katvealueiden syntyvan jopa yli 5 asteen mittauskulmilla. Keksinnon idean 
mukaisesti menetelmassa tulee kayttaS sellaista laserlaitetta, jolla voidaan lahettaa 
riittava maara laserpulsseja/naytteita yksittaisten puiden tai puuryhmien 
30 erottamiseksi toisistaan. Naytteiden maaraan vaikuttaa tietenkin haluttu tarkkuus, 
jolla puustoa halutaan mitata, ja puiden latvuston koko. Useimmissa tapauksissa 
naytteita kannattaa kerata useita neliometria kohden ja naytemaaran kasvattaminen 
kymmeneen tai jopa kymmeniin tuo selkea tarkkuusparannusta ainakin pohjoisessa 
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havumetsavyohykkeessa. Pulssi/naytemaSran optimointi on optimointia tarkkuuden 
ja kustannusten valilla. Riippuen edella mainituista parametreista on va.kea antaa 
ohjearvoa tarvittavalle pulssintoistotaajuudelle, mutta alarajana voidaan p.taa 50 
kHz yksittaiset puut nakyvat todella hyvin mikali vastaava arvo on no.n 200 kHz. 
Talis hetkella on olemassa vain muutama laite, jotka tayttaa namS vaatimukset. 
Esimerkiksi TopoSys-1 laserkeilaimen pulssintoistotaajuus on 80 kHz, mittaust.heys 
nimelliselta 800 m korkeudelta on 3-1 pulssia neliometriS kohden, ja mittaustarkkus 
x- ja y-suunnassa on parempi kuin 1 m ja tasaisilla a.ueiHa z-suuntainen tarkkuus on 
parempi kuin 15 cm. Tyypiliinen pulssitoistotaajuus muilla laitteilla on n. 10 kHz, 
mutta laitteistot kehittyvat koko ajan. 

Metsamittauksiin vaikuttaa myas vuodenaika. Talviaikaan (lunta maassa) 
iasermittauksia ai kannata suorittaa, koska lumi saattaa heijastaa pu.ssin po,spa,n 
,aserin ilmaisinosasta. Optimaalinen aika on silk*, kun latvustossa on tarpeeks, 
naulasia ja lehtia heijastaaksaen laserin signaalin latvasta ja mean pinnassa a, dm 
runsasta aluskasvillisuutta. Mittauksat Suoman o.osuhteissa ovat kuitenk,n 
osoittaneet, etta hyvalaatuista puustotietoa voidaan tuottaa myos 
toppukesasta/alkusyksysta. Radiotaajuisilla tutkilla mittauksia voidaan suonttaa 
kaikkina vuoden aikoina. Mikali taajuus on kuitenkin pienempi ku,n 5 GHz, 
lehdettemien lehtipuiden pituus tulee aliarvioitua. 

Koska mittauksat tulee tehda pienilla mittauskulmilla (vartikaaliin nahden), on 
mitattavan alueen leveys kapea (800 m lentokorkeudelta 282 metria levea olettaen 
10 astetta maksimaaliseksi poikkeutukseksi), minka vuoksi sovelluksesta n.ppuen 
voidaan joutua lentamaan useita rinnakkaisia linjoja koko alueen kattam.seksi. 

Keksinnon edullisissa suoritusmuodoissa eri mallien luomisessa kaytetSan hyvaksi 
eri laserianestelmissa kaytettavia erilaisia tallannustapoja. Tyypillisesti joka,sessa 
jarjestelmassa on seka ensimmaisen ja jalkimmaisen kaiun tallannustavat (engl. first 
, p 0 /sa mode ja /est pulse mode). Ensimmiisen kaiun tallennustapa ilmaisee laseria 
,ahimpana olavien kohteiden etaisyydan ja jalkimmaisen kaiun tallennustapa pyrkn 
nakemaan esimerkiksi metsaisessa kohteessa maanpinnan. Keksinnossa hava,tt„n, 
etta koska molemmat tallennustavat sisaltavat fetoa sake puuston latvaosista etta 
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maanpinnasta, saadaan paras tulos, mikali mittaustulokset yhdistetaan. Samalla 
voidaan kasvattaa kaytettavissa olevien naytteiden maaraa. Tulevaisuudessa voi olla 
my6s mahdollista, etta naiden tallennustapojen lisaksi tai sijasta tallennetaan koko 
profiili, kuten ei-keilaavilla jarjestelmilla on aikaisemmin tehty (J. Hyyppa, M. 
Hallikainen, 1993, "A helicopter-borne 8-channel ranging scatterometer for remote 
sensing, Part II: Forest Inventory", IEEE Transactions on Geoscience and Remote 
Sensing, 31:170-179. Tama kuitenkin edellyttaa huomattavaa datan talletus- 
kapasiteettia. 

Talla tavoin koko alue tulee katettua pienilla keiloilla, joille voidaan maarittaa 
koordinaatit. GPS-jarjestelmasta johtuen mittauspisteet ovat WGS-84-jarjestelmassa 
(lisatietoja kiijasta GPS-paikanmSSritys (Markku Poutanen 1998)) ja ne pitaa 
muuntaa kohteelle haluttuun koordinaattijarjestelmaan ja projektioon (kuten ykj ja 
kkj) kayttaen olemassa olevia konversiofunktiota. Kunkin maan kartoitus- ja 
maanmittausorganisaatiot ovat kykenevii toimittamaan tarvittavat korversiotiedot. 
Tuloksena syntyy maastopistetiedosto, joka kasittaa yksittaisia x-,y- ja z- 
koordinaatteja maaritettyna paikalliseen koordinaatti- ja projektiojarjelmaan. 
Koordinaatistojen ja projektioiden maarityksesta voi lukea teoksesta KSytSnnOn 
geodesia (Martti Tikka 1985). 

Kolmiulotteisen puuston pituusmallin laskeminen mitatusta aineistosta 

Pistetiedostoa analysoimalla voidaan tuottaa erilaisia pintamalleja. Seuraavassa on 
esitetty maaritelmat erilaisille malleille, koska useimmille niista ei ole virallista 
maaritelmaa. 

Pisteet, jotka ovat heijastuneesta maasta, synnyttavat maastomallin. Laserkeilaimen 
tuottamaa dataa ennen prosessointia kutsutaan yleensa pintamalliksi (DSM, digital 
surface model), josta voidaan sitten prosessoinnilla tuottaa maastomalli tai 
korkeusmalli (DTM, digital terrain model ja DEM, digital elevation model). 
Valitsemalla pisteet, jotka ovat heijastuneet kasvillisuuden korkeimmista kohdista, 
saadaan kasvillisuusmalli, jota voi kutsua myos metsaisilla alueilla latvamalliksi. 
Keksinnon mukaan on havaittu, etta mittaamalla puustosta etaisyystietoa riittavalla 
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maaralla mittauspisteita on mahdollista laskea pituusmalli, jossa on huomioitu 
maaston vaihtelevuus, ja josta siten voidaan johtaa puiden todellinen pituus toisin 
kuin latvamallissa, joissa ilmenee ainoastaan puiden korkeus tiettyyn 
referenssipintaan nahden. Pituusmalli voidaan esim. laskea latvamallin ja 
5 maastomallin/korkeusmallin erotuksesta, joka on siis metsaisilla alueilla 3-ulotteinen 
esitys siella esiintyvan puuston pituudesta. Mallit esitetaan joko hajapisteina tai 
saannollisena ruutuverkkona. 

Keksinto hyodyntaa juuri taman puuston pituusmallin kayttOa puuston inventoinnissa. 
10 Seuraavassa kuvataan esimerkki eraasta toimivasta ratkaisusta. 

Laskettaessa latvamallia, maastomallia tai puuston pituusmallia kannattaa toimia 
seuraavien periaatteiden mukaisesti. Muunnetaan mitatut pisteet halutulle ruudukolle 
(maaritetaan yksittaisen ruudun (jatkossa kutsutaan pikseliksi) koko ja alueen koko), 

15 kts. kuvio 3. Muunnettaessa pisteita ruudukkoon tulee muistaa, etta muunnos 
saattaa aiheuttaa virheita, varsinkin jos naytteita ei ole tarpeeksi kutakin ruutua 
kohden. Keksinnossa on havaittu, etta haluttaessa inyentoida puita yksitellen tai 
pienina ryhmina tulee pikselikoon olla esimerkiksi 1 m tai mielellaan 0.5 m, jolloin 
latvuksen leveyden maarittaminenluotettavasti on mahdollista. 0.5 m 

20 pikselikoko/erotuskyky vaatii naytteita runsas kymmenkunta neliOmetria kohden. 
Yksinkertainen tapa aloittaa maastomallin ja latvamallin laskeminen on valita 
kuhunkin pikseliin tulevista z-arvoista minimi ja maksimi. Saadaan kaksi pintaa, joista 
toinen kuvaa maastomallin maksimiarvoja (maksimipinta) ja toinen miminiarvoja 
{minimipinta). Maksimiarvot kuvaavat hyvin latvuston ylimpia osia ja aukeilla alueilla 

25 maaston pintaa. Minimiarvot kuvaavat maaston korkeusmallia, mutta joukossa on 
mm. puiden latvoja. Minimisuodattimen avulla voidaan minimiarvoista laskea nk. 
karkea maaston korkeusmalli. Suodattimen ikkunan koko tulee olla suurempi kuin 
yksittaisten latvojen halkaisija. KeksinnOssa havaittiin, etta hyva oletusarvo 
minimisuodattimen koolle on Suomessa 8 m runsaspuustoisilla alueilla. 

30 Minimisuodatuksen jalkeen tata nk. karkeaa maaston korkeusmallia verrataan 
alkuperaisien hajapisteiden koordinaatteihin. Pisteet, jotka poikkeavat tasta 
minimisuodatetusta pinnasta tietyn maaran D, esimerkiksi 1 metrin verran 
maksimissaan, luokitellaan maasta aiheutuneiksi pisteiksi. Naiden pisteiden avulla 
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interpoloidaan puuttuvien alueiden maaston pinta esimerkiksi Delaunayn kolmioinnin 
avulla. Syntynytta maaston korkeusmallia voidaan parantaa iteraation avulla 
seuraavasti: olettamalla tama pinta uudeksi referenssiksi ja verrattuna alkuperaisiin 
hajapisteisiin valitaan uudet maasta aiheutuneet pikselit, jolle operaatio toistetaah. 
5 Toistoa voidaan jatkaa, kunnes korvautuneiden pisteiden maara on vahamen ja laatu 
on tyydyttava. Toinen mahdollinen menetelma maastomallin maarittamiseen on 
esitetty julkaisussa (K. Kraus, N. Pfeiffere, 1998 "Determination of terrain models in 
wooded areas with airborne laser scanner data", ISPRS Journal of Photogrammetry 
& Remote Sensing, 53: 1993:203). 

10 

Mikali laserdatassa on reikia tai yksittaisia poikkeavia pisteita, voidaan ne korvata 
interpoloimalla sille uusi arvo viereisten pisteiden avulla. Yksittaisten poikkeavien 
pikselien tunnistamiseen voidaan laskea jokaiselle pikselille mm. gradientti. 
Poikkeavat arvot voidaan nain poistaa. 

15 

Kuviossa 4 on esitetty laserdatasta laskettu minimipinta, jossa kunkin pikselin sisalle 
kuuluvista pisteista on valittu z-arvoltaan pienin. Vastaava maksimipinta on esitetty 
kuviossa 5. Molempien kuvien pikselikoko on 50 cm ja kuva kattaa 125 m. Kuviossa 
6 on esitetty edella kuvatulla periaatteella (D = 1m), ilman iteraatiokierroksia ja 
20 delaunayn kolmioinnilla tehty maastomalli. Kuviossa 7 on esitetty vastaava puuston 
pituusmalli kyseisella alueella. 

Puustotunnusten maarittaminen puuston pituusmatlista 

25 Keksinnossa on havaittu ensimmaista kertaa maailmassa, etta hyvaresoluutioisesta 
puuston pituusmallista saadaan maaritettya seuraavat yksittaisen puun tunnukset 
ylhaalta nakyville puille: puun pituus, puun sijainti, latvuksen lapimitta, rungon 
lapimitta esim. 1.3 metrin korkeudelta, kehitysluokka, ika, runkotilavuus, pohjapinta- 
ala ja puulaji. Isommalle puurykelmalle voidaan lisaksi maarittaa naiden tunnusten ja 

30 naiden tilastollisen tiedon lisaksi ylhaalta nakyvan puuston runkoluku ja 
latvuspeittoprosentti, joita molempia voidaan kayttaa koealojen ja kuvioiden 
tunnuslukujen estimoinnissa. Aikaisemmissa tutkimuksissa ei ole pystytty 
hyOdyntamaan laserin avulla latvuspeittoprosenttiiin liittyvia suureita 
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profiilimittauksissa; myos talta osin keksinto tuo merkittSvaa uutta nykyiseen 
tietamykseen. Naiden tietojen avulla voidaan sitten maarittaa puuston om.na.suuden 
halutulta alueelta. 

Seuraavassa esitetaan esimerkki eraasta tavasta naiden tunnusten laskemiseksi. 

Laskennan ja estimoinnin lisaksi menetelma vol hy6dyntaa olemassa olevaa 
tietamysta, mika parantaa estimaattien luotettavuutta merkittavasti. Myos talta os.n 
keksintd on ainutlaatuinen. 

Puun rajaaman alueen sisalta saatujan naytepisteiden avulla maaritataan joko 
kaksiulottaisasti (poikkilalkkaus) tai kolmiulottaisasti puun gaomatria jafta, puun 
hahmo joko puulajin tunnistamiksi tai puuston mallintamiseksi. 

Puustotunnustan laskaminan puuston pituusmallista tapahtuu sauraavasti. 
Miaritataan yksittaistan puidan sijainti ja latvuston aariraunat parinteisilia hahmon- 
tunnistuksanmanetalmilia. Ylaisasti puidan sijaintian maaritys parustuu pa.kall.stan 
maksimian atsimisaan. Paikallisan maksimin ioytamisaksi lasarilla tuotattua 
pituuamallia suodatetaan alipaastasuodattimalla. Suodatuksan maaraan va,kuttaa 
kaytatty puuston pituusmallin pikselikoko ja aineiston lasernayttaidan lukumaara 
naliamatria kohdan. Tyypillisasti yksi suodatuskarta on MM "man suodatusta 
saattaa kuvalta IPytya samastakin latvuksesta useita paikallisia maks.maja. 
Latvuksan koon tai san aariraunojan maarittamiseen I6ytyy usaita asiaan sovaltuvia 
algoritmaja hahmontunnistuksan alalta. Latvuksian aarirajojan atsimisassa (ta, puun 
latvusta vastaavan aluaan maarittamisassa, sagmantoinnissa) paramatrit tulaa valrta 
siten atta kaksi arillista puuta ai tula yhdistattya yhdaksi latvaksi, mikaii va,n 
mahdollista. Paramatrian valinta sovaltumaan kaytossa olavalla puustotyyp,lla tulaa 
tend a manuaalisasti. Kuvio 8 on asimarkki au.omaat.isan 
hahmonttinnistusrnenetelman tuottamasta latvuksen rajoista. 

Latvuksan kaskimaarainan lavays L laskataan latvuksan paittaman pinta-alan A 
avulla. 



16 



112402 




Puun pituudeksi h oletetaan latvuksen sisalta loytyva puuston pituusmallin 
maksimipiste (puun korkein kohta). 

h=max(hi) (2) 

missa h, ovat yksittaisia puuston pituusmallin antamia pituuksia latvuksen pinta-alan 
A sisalla. Keksinnossa on selvitetty koemittauksin, etta yksittaisten puiden ja 
valtapuiden pituus voidaan maarittaa 1 m keskivirheella, kuvio 9. Vertailuaineistona 
oli hypsometrilla tehdyt maastomittaukset 89 puulle. Mittauksen harha oli 14 cm, 
mika on erittain hyva tulos. Selitysaste on 0.97. 

Puun sijainti maaritetaan maksimipistetta vastaavan x- ja y-koordinaattitiedon avulla. 

Pohjoisella havumetsaalueella puiden latvuksen keskimaaraisella leveydella L on 
selkea yhteys puun lapimittaan d 

d = oL + /1 (3) 

missa kertoimet a ja p kalibroidaan mieluiten paikallisen maastoaineiston avulla. 
Kertoimet tulee maarittaa kullekin puulajille (manty, kuusi, koivu, muu lehtipuu) 
erikseen, mikali jarjestelma hyQdyntaa puulajitiedon. 

Koska laserkeilain tuottaa myos pituustiedon ja koska lapimitalla ja pituudella on 
selkea yhteys, voidaan puulle mitattua pituutta kayttaa myos estimoinnissa. 
Kaytettaessa lineaarista regressiota puun lapimitan malli on muotoa. 

d = aL + J3h + r ( 4 > 
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missa kertoimet a, (5 Ja y kalibroidaan mieluiten paikallisen maastoaineiston avulla. 
Pituuden ja lapimitan suhde ei ole aivan lineaarinen, silla lapimitan kasvaessa 
riittavan suureksi kayran kulmakerroin pienenee lapimitan kasvaessa. Lisaksi my6s 
lapimitan kasvaessa pituuden hajonta kasvaa. Tarkempaa pituusmallia laadittaessa 
on siis ensin linearisoitava pituuden ja lapimitan valinen riippuvuus seka vakioitava 
pituuden varianssi lapimitan suhteen. Naslund (1936) (M. Naslund, 1936, 
-Skogsforseksanstaltens gallringsforsok i tallskog", Meddelanden fran Statens 
skogsfGrsOkanstalt 28 (1)kuvasi pituuden ja rinnaskorkeuslapimitan (lapimitta 1.3 m 
korkeudelta) riippuvuuden yhtalda seuraavasti 

A = T +1.3 (5) 

(a, +a 2 dy 

Kayttamalla epalineaarista muunnosta voidaan mallia (4) parantaa. Myos 
valtakunnallisten kayrien kayttaminen on mahdollista. Suomen osalta vastaavat 
kayrat ISytyvat viitteesta (J. Laasasenaho, 1982, 'Taper curve and volume functions 
for pine, spruce and birch", Communicatees Instituti Forestalls Fenniae 1 08, 74 s). 

Latvuksen lapimitan ja sita kautta yksittaisen puun lapimitan maarittamisen 
mahdollisia virhetekijoita on se, etta useita latvuksia tulee tulkittua kuuluvan samaan 
puuhun, jolloin latvuksen peittama pinta-ala ja saatu lapimitta ovat runsaasti oike.ta 
arvoja suurempia. Saatu latvuksen lapimitta L tuleekin tarkistaa esimerkiksi kaavan 
(3) avulla, josta saadaan karkea arvio rungon lapimitalle d. Mikali kaytossa on 
tilastotietoa siita, mika on todennakoisyys kullakin pituus- ja lapimittayhdistelmalla, 
voidaan poistaa useimmat virheellisesti maaritetyt latvuspinta-alat. Esimerkiksi mikal. 
esiintymistodennakSisyys on alle 1%, on todennakoista, etta puun latvuksen 
keskimaarainen leveys on maaritetty vaarin. Samaten olemassa olevaa tietamysta 
puun latvuksen keskileveydesta tulee kayttaa. Esimerkiksi puiden latvuston leveys 
on harvoin yli 8 m. Tallaiset tietamykseen perustuvat saannot eivat ole 
vaittamattomyys, mutta parantavat laskettujen puustotunnusten estimaatteja. Koska 
, laseraineistosta pystytaan maarittamaan seka yksittaisen puun pituus ja latvuksen 
keskilapimitan avulla rungon lapimitta, voidaan tietamykseen ja tilastotietoon 
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perustuvilla malleilla huomattavasti parantaa tilavuusestimaattien luotettavuutta, 
mika tekee keksirinOsta osaltaan ainutlaatuisen kaukokartoitusmenetelman 
maailmassa. lima- tai videokuviin perustuvat menetelmat pystyvat maarittamaan 
vain latvuston keskilapimitan. Ilman muuta mitattavaa suuretta vastaavia tarkistuksia 

5 ei voida toteuttaa nailla aineistoilla. Ja koska yksittaisen puun tilavuus maaritetaan 
perinteisesti tarkemmin pituuden ja lapimitan avulla kuin pelkan lapimitan avulla, on 
keilaavan laserin avulla mahdollisuus tuottaa huomattavasti tarkempia 
tilavuusestimaatteja kuin ilmakuvien tai videokuvien avulla. Alustavat tulokset 
viittaavat myos parempaan tarkkuuteen kuin perinteissa kuvioittaisessa 

10 inventoinnissa saavutetaan. 

Tietamysta voidaan kayttaa estimoimaan valtajakson alle jaavan puuston osuuden 
maaraa, silla laserkeilain tuottaa pituusmallia vain ylhaalta nahtavissa olevasta 
puustosta, Niinpa puustotunnukset saadaan maaritettya vain metsikossa esiintyvalle 
15 nakyvalle jaksolle. Mikali on olemassa tietamysta puuston todennakfiisesta 
pituusjakauma/lapimitta-jakaumasta, voidaan laserkeilain tuloksia parantaa 
kayttamalla naita jakaumia hyvaksi. Weibull-jakaumasta on tullut yleisimmin kaytetty 
funktio kuvaamaan lapimittajakaumaa. Suomessa Weibull-jakaumaa kaytettiin 
ensimmaisen kerran tutkimuksessa (P. Kilkki, R. Paivinen, 1986, "Weinbull-function 
20 in the estimation of the basal ansa DBH-distributions", Silva Fennica, 20: 149-156.). 
Weibull-jakauma tayttaa lapimittajakaumille tyypillisesti asetetut joustavuus- 
vaatimukset (A. Kangas, J. Kangas, K Korhonen, M. Maltama, R. PSivinen, 1990, 
Metsas kuvaavat mallit, Silva Carelica 17, 143 s.). Beta-jakauma on viela Weibull- 
jakaumaa joustavampi, mutta sen laskennallinen kasittely vaatii numerista 
25 integrointia skaalauskertoimen maarittamiseksi. Jakaumien kaytdssa verrataan 
laserin avulla maaritettya jakaumaan tiedossa olevaan puuston jakaumaan ja 
puuttuvien puiden osalta lasketaan korjaustekijat kullekin puustotunnukselle. 
Jakaumien kaytt6 ei ole valttamatdnta, mutta on suositeltavaa silloin, kun puuston 
tiedetaan olevan varsin tiheaa, todellinen jakauma on tiedossa ja pyritaan 
30 mahdollisimman pieneen systemaattiseen virheeseen, silla ilman korjausta 
laserkeilainaineisto tuottaa aliarviota puuston runkotilavuudesta. 
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Lapimitan ja pituuden avulla voidaan maarittaa puun kehitysluokka ja estimoida puun 
ika. 

Yksittaisen puun pohjapinta-ala g (m 2 /ha) saadaan kaavasta 

Runkoluku voidaan maarittaa yksinkertaisesti kuvasta maaritettyjen latvuksen 
lukumaarana. 

Latvuspeittoprosentti voidaan maarittaa latvuksien peittaman pinta-alan suhteena 
koko pinta-alaan. Latvuspeittoprosentti voidaan maarittaa myos siten, etta lasketaan 
puihin osuvien laserpulssien lukumaara suhteessa koko pulssimaaraan haiutulla 
alueella. 

Puulajin maarittamiseen on kaksi kayttOkelpoista menetelmaa. Julkaisussa J. 
Hyyppa, Hyyppa, H., Samberg, A., 1999, Assessing Forest Stand Attributes by Laser 
Scanner, Laser Radar Technology and Applications IV, Proceedings of SPIE, 3707: 
57-69 on esitetty laserkeilaimella tuotettu poikkileikkauskuva metsasta. Kartiomaiset 
kuusikot erottuvat helposti mannikoista ja koivikoista. Analysointi voidaan tehda 
myos 3-ulotteisesti, jolloin yhden puun rajaaman alueen sisalta saatujen 
naytepisteiden avulla maaritetaan puun geometria ja puun hahmo, jonka avulla se 
pyritaan tunnistamaan. Ilmakuvaus on verraten halpa media, joten ilmakuva 
yhdistettyna laserkeilainkuvaan antaa paremman tiedon puulajista. Talloin ilmakuva 
voidaan viela orto-oikaista laserilla maaritetyn puuston pituusmallin avulla (esim. 
suodatettu pituusmalli), mika on myos uutta maailmassa. TSIIoin saadaan puuston 
pituusmallissa ja ilmakuvassa olevat latvuksen saatua paallekkain. llmakuvassa 
(vaaravarikuva) lehtipuut erottuvat helposti mannikoista ja kuusikoista, joiden 
valillakin esiintyy selkea variero. Mikali laserlaite tuottaa intensiteettitietoa 
sironneesta tehosta, siita on kayttoa tehokkuutta puulajien luokitteluun. 
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Yksittaisen puun tilavuuden estimointi tapahtuu edella arvioitujen parametrien avulla. 
KSytettSvissa on kolme eri vaihtoehtoa: 1) tilavuuden estimointi pelkSn puuston 
pituuden avulla, 2) tilavuuden estimointi puuston pituuden ja maaritetyn lapimitan 
avulla ja 3) tilavuuden estimointi pituuden, lapimitan ja puulajin avulla. Laasasenaho 
5 (1982) (J. Laasasenaho, 1982, " Taper curve and volume functions for pine, spruce 
and birch", Communicationes Instituti Forestalls Fenniae 108, 74 s.) on esittSnyt 
funktiot, miten pituuden ja lapimitan avulla voidaan laskea yksittaisen puun tilavuus 
kullekin puulajille seuraavasti: 

10 mannylle v = 0.036089 d Z0139S (0.99676) d h 207025 (h -1 .3) - 107209 (7) 

kuuselle v = 0.022927 d 191505 (0.99146) d h 2 * 82541 (h -1.3) - 1 * 53547 (8) 

koivulle v= 0.01 1 1 97 d 2 - 10253 (0.98600) d h 3 * 98519 (h -1 .3) - Z6S900 (9) 

15 

Mallien keskivirheet ovat noin 7-8,5 %. 

Koealaa tai kuviota vastaavat puustotunnukset saadaan maaritettya laskemalla 
vastaavat tunnukset kSyttaen kaikkia yksittaisen puun tunnuksia halutulla alueella. 

20 Esimerkiksi kuviota vastaa tilavuus per hehtaari voidaan laskea summaamalla 
kaikkien yksittaisten puiden tilavuuden kyseisen kuvion sisalta ja jakamalla tilavuus 
kuvion koolla. Aikaisemmin esitetylla jakaumaperiaattella voidaan parantaa 
estimaattien tarkkuutta nakymattomSn puuston osan osalta. Talloin tosin tulee 
muistaa, etta puuston tilastollisesta jakaumasta on oltava tietoa. Keskipituus 

25 saadaan joko aritmeettisena keskiarvona yksittaisten puiden pituuksille. YleensS 
metsien inventoinnissa on tapana laskea painotettu keskipituus, jossa 
painokertoimena on kunkin yksittaisen puun pohjapinta-ala. Pohjapinta-ala hehtaaria 
kohden saadaan summaamalla kaikkien puiden pohjapinta-alat ja jakamalla se 
kuvion koolla. 

30 

Menetelman yksinkertaistettu tapa tuottaa tilavuusinformaatio on arvioida latvuksien 
koon penjsteella latvuspeittoprosentti, pohjapinta-ala tai vastaava tiheyssuure ja 
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puuston keskipituuden (maaritetty puuston keskipituusmallin avulla) maarittaa 
puuston tilavuus llvessalon relaskooppitaulukkojen tapaan (Tapion Taskukirja 1978). 

Useimmat metsat on jo useampaan kertaa inventoituja. Esitetyn menetelman 
tuloksia parantaa vanhan inventointitiedon kaytt6 olettaen, etta vanha inventointitieto 
on luotettavaa. Vanhaa inventointietoa tulee kayttaa siten, etta valtetaan virheet, 
jotka ovat todellisuudessa mahdottomia. Esimerkiksi kuvio, joka 10 vuotta sitten oli 
tilavuudeltaan 5 m^ha ei voi olla nyt tilavuudeltaan 200 m 3 /ha. 

Tehdyt koetyot ovat osoittaneet, etta kehitetty menetelma on tarkkuudeltaan 
vertailukelpoinen ja jopa parempi kuin perinteisesti kaytetyt operatiiviset inventointi- 
menetelmat. 

Menetelma soveltuu yksittaisten puiden, puuryppaiden, koealojen, kuvioiden ja 
laajempien alueiden inventointiin. Seuraavassa on esitetty menetelman 
soveltuvuutta kuvioittaiseen arviointiin ja valtakunnan metsien arviointiin, jotka ovat 
kaksi taloudellisesti tarkeinta inventointitapaa Suomessa. 

Kuvioittaiseen arviointiin menetelma soveltuu erinomaisesti erityisesti, jos se 
yhdistetaan vanhaan jo kerattyyn kuvioaineistoon. Talldin laserin avulla lasketut 
puustotunnukset voidaan yhdistaa muihin alueelta tarvittaviin tietoihin, kuten 
kasvupaikkatietoihin, jotka eivat ajan kanssa muutu kovinkaan paljoa. Laserin avulla 
maaritetya maastomallia voidaan kayttaa my6s avainbiotooppialueiden 
kartoittamisessa. 

Valtakunnan metsien inventointi on esitetylle menetelmalle erittain sovelias 
sovelluskohde. Laserkeilaimella yleensa tuotetaan kohteesta muutama sata metria 
levea kaistale. Haluttaessa peittaa laaja alue joudutaan lentamaan useita 
rinnakkaisia linjoja. Linjojen vaiisiin kaarroksiin, joissa differentiaali-GPS ei saa 
kadottaa satelliitteja, kuluu kuitenkin varsin paljon aikaa, mika nostaa menetelman 
kustannuksia. Valtakunnan metsien inventoinnissa laserkeilaimen kaltainen aineisto 
mahdollistaa sen, etta maastotyon maara vahenee ja laserkeilaimen avulla laskettuja 
yksittaisten puiden ja koealojen tietoja voidaan kayttaa satelliittiaineiston 
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opettamiseen. Valtakunnan metsien inventoinnissa myos lentaminen on tehokasta. 
Linjat voidaan suunnitella kattavan koko valtakunnan. Menetelman avulla 
valtakunnan jnetsien inventointi on mahdollista automatisoida miltein kokonaan. 
Metsantutkimuslaitoksessa on pyritty kehittamaSn kuvaavan spektrometrin (Aisa) 
avulla puustotiedon keraamista ja koealatiedon vShentamista. Aisalla saadut tulokset 
ovat kuitenkin merkittavasti huonommat kuin kayttamalia taman keksinnon 
menetelmia. 
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PATENTTIVAATIMUKSET 

1. Menetelma puustotunnusten m§arittamiseksi puuston yiapuolella mittaavan 
laitteen avulla, 

tunnettu siita, etta 

a) puustosta kerataan laserkeilaifnella kolmiulotteista tietoa kayttamalla 
sellaistamaaraa naytepisteita, etta yksittaiset puut tai puuryhmat erottuvat 
toisistaan, 

b) laserkeilaimella ker§tyista tiedoista tuotetaan kolmiulotteinen pituusmalli 
puuston pituudesta, 

c) saadusta pituusmallista maaritetaan puustotunnuksia, jotka ovat yksittaisten 
puiden tai puuryhmien ominaisuuksia ja/tai niiden avulla saatuja laajempien 
aiueiden puustotietoja. 

2. Patenttivaatimuksen 1 mukainen menetelma, tunnettu siita, etta 
naytepisteiden maara (pulssintoistotaajuus) valitaan halutun tarkkuuden mukaan. 

3. Patenttivaatimuksen 1 tai 2 mukainen menetelma, tunnettu siita, etta 
kohdassa a) keratyt tiedot ovat kolmiulotteisia etaisyysmittauksiin perustuvia 
tietoja. 

4. Jonkin patenttivaatimuksen 1-3 mukainen menetelma, tunnettu siita, etta 
kolmiulotteinen esitys puuston pituudesta (puuston pituusmalli) saadaan 
laskemalla pintamalleista puuston pituusmalli puuston ylimpia osia vastaavan 
latvamallin seka maanpintaa vastaavan maastomallin erotuksesta. 

5. Jonkin patenttivaatimuksen 1 - 4 mukainen menetelma, tunnettu siita, etta 
kolmiulotteisen puuston pituusmallin tuottamisessa kaytetaan hyvaksi erilaisia 
tallennustapoja ja yhdistamalla naiden mittaustuloksia, jolloin kaytettavissa 
olevien naytepisteiden maara kasvaa. 

6. Jonkin patenttivaatimuksen 1 - 5 mukainen menetelma, tunnettu siita, etta 
puun rajaaman alueen sisaita saatujen naytepisteiden avulla maaritetaan joko 
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kaksiulotteisesti (poikkileikkaus) tai kolmiulotteisesti puun geometria ja/tai puun 
hahmo joko puulajin tunnistamiseksi tai puuston mallintamiseksi. 

7. Jonkin patenttivaatimuksen 1-5 mukainen menetelma, tunnettu siita, etta 
5 yksittaisten puiden tai puuryhmien ominaisuudet, jotka saadaan pituusmallia 
analysoimalla ovat puiden sijainti, ika, pituus, latvuksen lapimitta, rungon 
lapimitta, pohjapinta-ala, kehitysluokka, puulaji, runkotilavuus ja/tai runkoluku 
pinta-alayksikkoa kohden ja naiden avulla johdettavat tilastolliset tunnukset. 

10 8. Patenttivaatimuksen 7 mukainen menetelma, tunnettu siita, etta yksittaisten 
puiden sijainti ja latvuksen aarireunat analysoidaan sinansa tunnetulla 
hahmontunnistusmenetelmallS. 

9. Patenttivaatimuksen 7 mukainen menetelma, tunnettu siita, etta puun pituus 
15 saadaan maarittamalla sinansa tunnetulla menetelmalla latvuksen sisalta loytyva 

puuston pituusmallin maksimipiste. 

10. Patenttivaatimuksen 7 mukainen menetelma, tunnettu siita, etta puun rungon 
lapimitta voidaan johtaa latvuksen keskilapimitan tai puun pituuden ja latvuksen 

20 keskilapimitan avulla ja mahdollisesti kayttamalla hyvaksi tietamykseen 
peoistuvia saantoja ja mahdollisesti kullekin puulajille erikseen. 

1 1 . Patenttivaatimuksen 7 mukainen menetelma, tunnettu siit§, etta puun 
kehitysluokka ja ika maaritetaan lapimitan ja pituuden avulla, ja mahdollisesti 

25 kullekin puulajille erikseen. 

12. Patenttivaatimuksen 7 mukainen menetelma, tunnettu siita, etta runkoluku 
maaritetaan kuvasta maaritettyjen latvuksien lukumaarana. 

30 13. Jonkin patenttivaatimuksen 1-5 mukainen menetelma, tunnettu siita, etta 
latvuspeittoprosentti maaritetaan latvuksien peittaman pinta-alan suhteena koko 
pinta-alaan. 
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14. Jonkin patenttivaatimuksen 1-5 mukainen menetelma, tunnettu siita, etta 
latvuspeittoprosentti maaritetaan s'rten, etta lasketaan puihin osuvien 
laserpulssien lukumaara suhteessa koko pulssimaaraan halutulla alueella. 

5 15. Jonkin patenttivaatimuksen 8-14 mukainen menetelma, tunnettu siita, etta 
isommalle puurykelmalle maaritetaan yksittaisten puiden tai puuryhmien 
ominaisuuksien ja naiden tilastoilisen tiedon lisaksi ylhaalta nakyvan puuston 
runkoluku ja latvuspeittoprosentti, joita molempia voidaan kdyttaa koealojen ja 
kuvioiden tunnuslukujen estimoinnissa. 

10 

16. Jonkin patenttivaatimuksen 8-15 mukainen menetelma, tunnettu siita, etta 
puuston tilavuus maaritetaan kokonaan tai osittain puuston keskipituuden ja 
latvuspeittoprosentin (latvuksien osuus kokonaispinta-alasta)/pohjapinta-alan 
avulla. 

15 

17. Jonkin patenttivaatimuksen 1-16 mukainen menetelma, tunnettu siita, etta 
puustotunnusten maarittamiseen kaytettyihin tietoihin yhdistetaan muista 
menetelmista saatuja tietoja, kuten vanhaa inventointitietoa alueelta, ilmakuvia ja 
satelliittikuvia, laser-keilaimella saatujen tietojen lisaksi. 

20 

18. Jonkin patenttivaatimuksen 1-17 mukainen menetelma, tunnettu siita, etta 
kolmiulotteisesta puustotiedosta maaritetaan puun laatuun liittyvia tunnuksia 
maarittamaiia oksattoman rungon pituus. 

25 19. Jonkin patenttivaatimuksen 1-18 mukainen menetelma, tunnettu siita, etta 
puustotunnusten maarittaminen laserkeilaimella saadusta tiedosta tapahtuu 
tietokoneohjelman avulla. 

20.Tietokoneohjelma puustotunnusten maarittamiseksi tiedoista, jotka on saatu 
30 puuston yiapuoleila mittaavan laitteen avulla, 

tunnettu siita, etta se suorittaa seuraavat vaiheet 
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a) kolmiulotteisen puuston pituusmallin tuottaminen puustosta keratysta 
kolmiuiotteisesta tiedosta puustoalueelta, jossa on kaytetty sellaista maaraa 
naytepisteita, etta yksittaiset puut tai puuryhmat erottuvat toisistaan, 

b) puustotunnusten maarittaminen saadusta pituusmallista, jotka 
5 puustotunnukset ovat yksittaisten puiden tai puuryhmien ominaisuuksia ja/tai 

niiden avulla saatuja laajempien alueiden tietoja. 

21 . Patenttivaatimuksen 20 mukainen tietokoneohjelma, t u n n e 1 1 u siita, etta se 
suorittaa joinkin patenttivaatimuksen 4-18 mukaisen menetelmSn. 



10 



— 27 



112402 



PATENTKRAV 

1. Forfarande f6r att definiera kannetecken pa trddbestand med hjaip av en 
anordning som mater ovanom tradbestandet, 
5 kannetecknat darav att, 

a) tredimensionell information om tradbestandet samlas medelst en 
Iaseravs6kare genom andvandning av en sadan mangd provpunkter att 
enskilda trad eller tradgrupper kan sarskiljas, 

b) av data som samiats med Jbjajp av laseravsdkaren produceras en 
10 tredimensionell langdmodell av langden pa tradbestandet, 

c) av den erhallna langdmodellen defmieras kannetecken pa tradbestandet, vilka 
ar egenskaper pa individuella trad eller tradgrupper och/eller information om 
tradbestand Sver mera omfattande omraden. 

15 2. Fdrfarande enligt patentkrav 1, kannetecknat darav att mangden 
provpunkter (impulsfrekvensen) valjs efter 6nskad noggrannhet. 

3. Forfarande enligt patentkrav 1 eller 2.kannetecknat darav att data som 
samiats i steg a) baserar sig pa tredimensionella avstandsmatningar. 

20 

4. Fdrfarande enligt nagot av patentkraven 1-3, kannetecknat darav att 
man erhaller en tredimensionell framstailning over langden pa tradbestandet (en 
langdmodell av tradbestandet) genom att fran ytmodeller nlkna ut langdmodellen 
pa tradbestandet fran skillnaden mellan toppmodellen som svarar mot de Sversta 

25 delama av tradbestandet och terrangmodellen som svarar mot markytan. 

5. F6rfarande enligt nagot av patentkraven 1-4, kannetecknat darav att 
man vid produktionen av den tredimensionella langdmodellen av tradbestandet 
utnyttjar olika lagringsmetoder och kombinerar matresultat som erhSllits medelst 

30 dessa, varvid mangden provpunkter som star till anvandning 6kar. 

6. F6rfarande enligt nagot av patentkraven 1-5, kannetecknat darav att 
man med hjaip av provpunkter som erhallits innanf6r ett omrade som begransas 
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av tradet definierar tradets geometri antingen tvadimensionellt (tvarsnittet) eller 
tredimensionellt och/eller gestalten pa tradet for att antingen definiera tradslaget 
eller gora en modell av tradbestandet. 

5 7. Fdrfarande enligt nigot av patentkraven 1-5, kannetecknat darav att 
egenskapema pa enskilda trad eller tradgrupper som erhalles genom analysering 
av langdmodellen ar laget, ildem, langden, tvarsnittet pi kronan, tvarsnittet pa 
stommen, bottenarean, utvecklingsklassen, tradslaget, stomvolymen och/eller 
stomantalet per ytenhet och statistiska kannetecken som harleds fran dessa. 

10 

8. Fdrfarande enligt patentkrav 7, kannetecknat darav att laget pa 
individuella trad och ytterkanterna pa kronan analyseras med en i och f6r sig 
kSnd gestalt- identifieringsmetod. 

15 9. F6rfarande enligt patentkrav 7, kannetecknat darav att langden pa tradet 
erhalles genom att med en i och f6r sig kand metod definiera maximipunkten av 
langdmodellen av tradbestindet, vilken man kan finna innanfOr kronan. 

10. Ffirfarande enligt patentkrav 7, kannetecknat darav att tvarsnittet pa 
20 tradstommen kan harledas med hjalp av det genomsnittliga tvarsnittet pa kronan 
och/eller tradets langd och det genomsnittliga tvarsnittet pa kronan och eventuellt 
genom att utnyttja regler som baserar sig pa kunskap och eventuellt skilt fdr vart 
och ett tradslag. 

25 1 1 . Fdrfarande enligt patentkrav 7, kannetecknat darav att tradets 
utvecklingsklass och alder definieras med hjalp av tvarsnittet och langden och 
eventuellt skilt fdr vart och ett tradslag. 

12. Fdrfarande enligt patentkrav 7, k a n n e t e c k n a t dar av att stomantalet 
30 definieras som antal toppar som definierats fran bilden. 



29 



112402 



13. FOrfarande enligt nagot av patentkraven 1-5, kSnnetecknat darav att 
krontackningsprocenten definieras som fdrhallandet mellan arean som tacks av 
kronorna och ytan. 

14. FOrfarande enligt nagot av patentkraven 1-5, kannetecknat darav att 
krontackningsprocenten definieras pi sadant satt att antalet laseriulser som 
traffar traden i fOrhallande till hela pulsmangden raknas inom ett Onskat omrade. 

15. FOrfarande enligt nagot av patentkraven 8-14, kSnnetecknat darav att 
f6r en stdrre tradgrupp definieras utom egenskaper pa individuella trad eller 
tradgrupper och statistisk information dessutom stomantalet och krontacknings- 
procenten av tradbestandet som syns ovanrfran, vilka bSgge kan anvandas vid 
estimering av kannetecken pa provareor och figurer. 

16. FOrfarande enligt nagot av patentkraven 8-15, kannetecknat darav att 
volymen pa tradbestandet definieras helt eller delvis med hjalp av den 
genomsnittliga langden och krontackningsprocenten (andelen toppar av den 
totala arean)/ bottenarean. 

17. F6rfarande enligt nagot av patentkraven 1-16, kannetecknat darav att 
man till data som anvants fdr att definiera kannetecken pa tradbestindet fdrenar 
information som erhallits med andra metoder, sasom gammal 
inventeringskunskap fran omradet, luftfotografier och satellitbilder som laggs till 
data som erhallits med laseravsOkaren. 

18. FOrfarande enligt nagot av patentkraven 1-17, kannetecknat darav att 
man fran den tredimensionella informationen om tradbestandet definierar 
kannetecken i anslutning till tradslaget genom definiering av langden pa den 
grenfria stommen. 

19. FOrfarande enligt nagot av patentkraven 1-18, kannetecknat darav att 
definitionen av kannetecknena pa tradbestandet fran data som erhallits med 
laseranordningen sker med hjalp av ett datoiprogram. 
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20. Datorprogram f6r att definiera kannetecken pa tradbestand fran data som 
erhaliits med hjalp av en anordning som m§ter ovanom tradbestandet, k a n n e t 
e c k n a t darav att det utfSr fdljande steg 

a) producerar en tredimensionell langdmodell av tradbestandet fran 
5 tredimensionell information som samlats fran omradet av tradbestandet, dar 

man anvant ett sadant antal provpunkter att enskilda trad eller tradgrupper 
sarskils, 

b) definierar kSnnetecken pa tradbestandet fran den emallna langmodellen, 
vilka kannetecken pa tradbestandet ar egenskaper pa enskilda trad eller 

10 tradgrupper och/eller information om vidare omraden som ertiallits med hjSIp 

av dessa. 

21. Datorprogram enligt patentkrav 20, kannetecknat darav att det utf6r ett 
fSrfarande enligt nagot av patentkraven 4 

15 
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